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Within finite temperature field theory, we show that truncated non-perturbative self- 
consistent Dyson resummation schemes can be renormalized with local counter terms defined 
at the vacuum level. The requirements are that the underlying theory is renormalizable and 
that the self-consistent scheme follows Baym's ^-derivable concept. The scheme generates 
both, the renormalized self-consistent equations of motion and the closed equations for the 

P^ ' infinite set of counter terms. At the same time the corresponding 2PI-generating functional 

and the thermodynamical potential can be renormalized, in consistency with the equations of 
motion. This guarantees the standard <i>-derivable properties like thermodynamic consistency 
and exact conservation laws also for the renormalized approximation schemes to hold. The 
proof uses the techniques of BPHZ-renormalization to cope with the explicit and the hidden 

0^ , overlapping vacuum divergences. 

PACS numbers: ll.lO.-z, ll.lO.Gli, ll.lO.Wx 
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In recent years the question how to appropriately treat particles in a hot and dense medium 
has continuously gained growing interest in many areas of physics ranging from plasma physics, 
condensed matter physics to nuclear and particles physics. Within a Green's function formal- 
in ' ism a consistent treatment of such phenomena frequently leads to consider dressed propagators, 
which follow from non-perturbative Dyson resummation schemes, rather than perturbative ones, 






Qh. in particular, if damping width effects play a significant role. 

Already in the early sixties, based on a functional formulation of Luttinger and Ward Q and Lee 



and Yang [y], Kadanoff and Baym [^ considered a class of self-consistent Dyson approximations. 
Baym reformulated this in terms of a variational principle, defining the so called ^-derivable 
rS \ approximations Q . Since in principle the truncation of the diagrammatic series of the functional 

C^ ' $ can be at arbitrary level it gives rise to a variety of approximations including Hartree and 

Hartree-Fock as the simplest schemes. The main virtue of this concept is that the resulting 
equations of motion are conserving and the corresponding equilibrium limit is thermodynamically 
consistent. This functional treatment constitutes the basis for the two-particle irreducible (2PI) 
diagram technique, where the functional $ generates the driving terms for the equations of motion, 
like the self-energy. Later the concept was extended to the relativistic case and formulated within 
the path integral approach by Cornwall, Jackiw and Tomboulis [^. There is no formal problem 
to extend the formalism to the Schwinger Keldysh real-time path method |^, applicable to the 
general case of non-equilibrium many-body theory. 

Despite such early conceptual formulations most applications of self-consistent approximations 
were pursued on the Hartree or Hartree-Fock level sometimes supplemented by RPA resummations 
(see, e.g., [^, ^) or perturbative estimates of higher order corrections. Thus, essentially mean field 
corrections to the self-energies were considered. Genuine two-point or even multi-point contribu- 
tions to the self-consistent self-energy, which give rise to a finite damping width, imply a new level 
of complexity. Various new conceptual problems, like leaving the quasi-particle picture or the is- 
sue of renormalization come in with considerable complications for the numerical solutions of such 
problems. In the pioneering work of Bielajew and Serot [ll^] for the first time the renormalization of 



self-consistent two-point self-energy loops were investigated at zero temperature but finite matter 
density. 

In recent years with the special interest in dense hadronic matter problems the <l>-derivable 
schemes with higher order self-energy terms were used to derive transport equations |11, 1^, 13, 14] 



from the corresponding Kadanoff-Baym equations |15] for the consistent and also conserving |16] 
transport treatment of particles with finite spectral width beyond the quasi-particle approximation. 
Also first investigations involving finite mass width effects on vector mesons were investigated 
within a self-consistent Dyson resummation scheme |17]. In most of these cases, however, the 
question of renormalization was circumvented by taking into account the imaginary part of the 
self-energy only, while the real part was neglected, or cut-off recipes - mostly symmetry violating 
- were employed, or even the counter terms were chosen temperature dependent! Yet, especially in 
the study of phase transitions, e.g., within chiral hadronic models of QCD, or for non-perturbative 
corrections of hard thermal loop approaches to QCD, e.g., within a ^-derivable scheme [ p!8| , [l9| ] it 
is important to consistently take into account both, real and imaginary parts of the self-energies. 

Therefore in this paper we address the more formal question of renormalizability of such non- 
perturbative approximations. We essentially concentrate on the thermodynamic equilibrium case 
and show how to obtain finite self-consistent dynamical quantities like the in-medium equations of 
motion and the self-energy of the particles and thermodynamic quantities like the pressure and the 
entropy. For definiteness and clarity of the presentation we use the (p^-theoicy as the most simple 
example to study the related questions. The results and techniques can easily be transfered to 
other theories. 

The paper is organized as follows: In section II we briefly summarize Baym's ^-functional using 
the combined real and imaginary time contour appropriate for thermal equilibrium within the path 
integral formalism |^. 



In section III we derive the general formalism for the renormalization of the self-consistent 



self-energy at finite temperature and the in-matter generating functional T. With the help of 



Weinberg's convergence theorem |2C] and the BPHZ-formalism of renormalization theory |21, |2 



we show that, in close analogy to perturbative renormalization (see, e.g., ||2^, ^, 25|) any <l>- 
derivable self- consistent approximation scheme can be rendered finite by subtracting pure vacuum 
counter terms given by closed recursive equations. Indeed the main complication arises from the 
fact that the self-consistent propagator is involved in divergent loops which gives rise to "hidden" 
divergences which have to be resolved. This leads to a Bethe-Salpeter equation for the divergent 
vacuum pieces with a kernel compatible with the functional $, which needs to be renormalized. 
The renormalized equations of motion for the self-consistent propagator are shown to be consistent 
with the renormalized 2PI generating functional which proves the consistency of counter terms at 
both levels. The diagrammatical interpretation shows that in strict analogy to perturbative renor- 
malization of thermal quantum field theory this procedure can be interpreted as renormalization 
of the wave functions, the mass and the coupling constants in the vacuum. We also give a closed 
expression for the renormalized self-consistent thermodynamical potential. 

Numerical solutions for the renormalized self-consistent Dyson equations beyond the standard 
Hartree approximation up to the self-consistent sunset-diagram level could be achieved; the results 
are discussed in the second paper of this series 



II. ^-DERIVABLE APPROXIMATIONS 

In the case of thermal equilibrium the real and imaginary time formalism can be combined by 
extending the Schwinger-Keldysh contour Cr, running from ti to tf and back to tj, by appending a 
vertical part Cxh running from ij to — i/3 (see fig. Q). One uses the fact that the factor exp(— /3H) 
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FIG. 1: The Schwinger-Keldysh real-time contour modified for the application to thermal equilibrium of 
quantum field theory. 



in the canonical density operator can be formally treated as a time evolution in imaginary time 
direction. The functional integral for bosonic fields has to be taken over all fields fulfilling the 
periodic boundary condition (j){ti — \(3) = <j){ti) which leads to the Kubo-Martin-Schwinger (KMS) 



condition for the Green's functions |27|. Since the equilibrium state is invariant in time one can take 
ti — > — oo and tf -^ +00 which is convenient when formulating the theory in energy-momentum 
space via a Fourier transformation. 

In addition to the usually introduced one-point auxiliary external source also a two-point aux- 
iliary external source is included. Variations of the latter generate contour-ordered expectation 
values of the form {Tc4>{xi)(f)[x2)) where C denotes the extended Schwinger-Keldysh time contour, 
and Tc stands for the ordering of the operators due to the ordering of the time arguments along 
this contour. 

The corresponding generating functional is defined within the path integral formalism of quan- 
tum field theory as 



Z[J,B] =N B(l) exp 



iS[</.] +i / d(l) JicPi + - d(12) ^120102 



(1) 



where A'^ is a normalization constant chosen below. Here and in the following the shorthand 
notation 



/ d(12 . . . n) /i2...n = / d^xi • • • d'^xnfixi, . 

Jc JCxRd-^ 



(2) 



for integrations in d-dimensional space is used (in the sense of dimensional regularization^ ) . The 
time integration has to be performed along the time contour introduced above. It is clear that also 
the action functional S has to be defined as the C x K*^"^ -integral of the Lagrangian. 
The generating functional W for connected diagrams reads 



W[J,B] = -i\n{NZ[J,B]). 



(3) 



^ We use this notation for convenient regularization only, in order to write down sensible non-renormalized func- 
tionals. 



The mean field and the connected Green's functions are defined by 

The last formula immediately follows from the definition of the partition sum (|l|) and (|3|) using the 
Feynman-Kac formula 

(7^0(xi) • • • <t){xn)) = / D(/) (/)(xi) • • • (l){xn) exp(iS' [(/)]). (5) 

Since the real-time part of this contour is closed by itself it can be shown that the functional ([l|) 
factorizes into the real-time part and the imaginary time part. Thus the Feynman rules for the 
calculation of connected Green's functions apply separately to the real and the imaginary part of 
the contour |£^]) since the functional (y) splits into a sum of the contributions from the vertical 
and the real-time part of the contour respectively: 

Z[J,B] = Z4J,B]ZTh[J,B],W[J,B] = W4J,B] + VFxh [ J, S] - i In A^ . (6) 

Within the real-time part the effect of the heat bath is completely taken into account by means of 
the analytical properties of the Green's function. It is uniquely determined by the KMS-condition 
which itself is a consequence of the above mentioned periodic boundary conditions for bosons within 
the path integral (see appendix |^ for details about the analytic properties of Green's functions 
and self-energies). 

By a functional Legendre transformation in ip and G one obtains the effective quantum action: 

r[ip, G] = W[J, B]- J d(l) ipiJi -^f d(12) {ip^ip2 + iGi2)Si2 - iln iV. (7) 

Now, as is well known from the usual functional formalism of quantum field theory, a formal saddle 
point expansion of the effective quantum action is an expansion in orders of h around the classical 
solution, where G is considered as an independent quantity, gives 

r[^, G] = S[^] + ^ TVlnCM^G-i) + ^ / d(12) D^^\Gi2 - Du) + $[<P, G]. (8) 

Herein the free propagator in the presence of a mean field is given by 

^^ 6ipi6ip2 ' 

The arbitrary constant M^ account for the overall normalization and cancels for any physical 
quantities as we shall see below. 

In the case of an ideal gas it is sufficient to subtract the pure vacuum part to render this 
functional finite, which leads to the well known result. At T = this subtraction corresponds to 
the renormalization of the total ground-state energy to zero. 

As we shall discuss below in the case of interacting particles this description is not sufficient 
to render the effective action and thus the pressure finite, since we need additional subtractions of 
vacuum sub-divergences to renormalize it. 

In the above sense the functional $ in (^) contains the parts of order 0{h'^) and higher. Since 
simple power counting shows the h-oicder of diagrams to be identical with the number of loops, ^ as 
a functional of G consists of all closed diagrams with at least 2 loops. The lines within the diagrams 
stand for dressed Green's functions G while the vertices are the bare vertices of the classical action 



with presence of the background field (p which can be immediately read off from Si [(p + (j)'] around 
(/)' = beginning at order (j)'^. 

The equations of motion are now given by the fact that we like to study the theory with vanishing 
auxiliary sources j and B. From (^) and ^ we obtain 

<5r I (5r i I 

Using (^) the equations of motion read 

The first line is of the form of a Klein-Gordon equation with the quantum corrections to the 
classical theory on the right hand side. The second equation is the Dyson equation and shows that 
the variation of the ^-functional with respect to G is the self-energy defined with respect to the 
classical field dependent propagator (P). This shows that the ^-functional must be two-particle 
irreducible. No propagator line must contain a self-energy insertion. In other words the closed 
diagrams representing contributions to $ must not split in disconnected pieces when cutting two 
lines. Diagrammatically the derivative of a functional with respect to G corresponds to opening 
one line of the diagrams representing it. In that sense ^ is the generating functional for skeleton 
diagrams for the self-energy where the lines represent fully dressed propagators. Thus the functional 
formalism avoids double counting in a natural way by omitting all non-skeleton diagrams from the 
Dyson-resummed equations of motion. Altogether the ^-functional formalism provides a closed 
system of exact equations of motion for the full 2-point function and the full mean field. Solving 
these equations would be equivalent to finding the full propagator of the quantum field theory 
which of course is impossible in practice. 

One obtains approximations by truncating the series for <^ at a certain vertex or loop order 
(which corresponds to the respective order in the coupling X oi h respectively), while preserving 
the forms (110) and ( pT] ) of the self-consistent equations of motion. Approximations of this kind 
respect the conservation laws for the expectation values of Noether currents for symmetries which 
are linearly operating on the field operators (including space-time symmetries and the according 



conserved quantities as energy, momentum and angular momentum) y, 11]. 

In the case of thermal equilibrium, setting the mean field and the propagator to the so- 
lution of the self-consistent equations the effective action gives the grand canonical potential 
$7 = —Tr[ip,G]j^K=o = —TlnZ{(3) p7[| . Since the real-time part of the contour in figure ^ is 
closed it vanishes for the solution of the equations of motion. A short summary about the analytic 
continuation from the real to the imaginary time formalism is given in appendix ^. 

Thus the formalism leads to a well defined treatment for bulk thermodynamical quantities of 
the system (like energy, pressure, entropy, etc.). 

All these quantities can be calculated either with real-time Green's functions or with the corre- 
sponding imaginary time functions, because as summarized in appendix |A| real and imaginary-time 
propagators are connected by the analytic properties of the Green's functions originating from the 
KMS-condition. For our purpose the real-time formalism is preferred, because of its simplicity 
with respect to the analytic structure of Green's functions, which easily permits to deal with the 
mixture of finite temperature and vacuum pieces occurring in the subtraction scheme. This also 
avoids the necessity to perform an analytic continuation from imaginary time to real-time Green's 
functions which is complicated to obtain for numerical results. 



In order to exemplify the method we apply the formalism to (/i^-theory with the Lagrangian 

^ = l{d,<p){d^c^)-^c^'-^cp\ (12) 

III. RENORMALIZABILITY OF ^-DERIVABLE APPROXIMATIONS 

In this section we show that in close analogy to the renormalization of perturbative diagrams also 
any ^-derivable self-consistent approximation can be renormalized with help of local temperature- 
independent counter terms. 

The proof uses the same line of arguments as in the perturbative case: The reason is, that 



the renormalization theory completely rests on Weinberg's power counting theorem |2C] which is 
formulated for a general class of Green's functions with a given asymptotic behavior. It does not 
depend on the special form of the propagators. 

The first step is a simple topological argument leading to the superficial degree of divergence 
for a given diagram 7, which for the (p^ theory simply follows from the number E of external lines 



ii 



5(7) =4-E. (13) 

Due to field reflection symmetry only diagrams with an even number of external lines are different 
from 0, the only divergent functions are those represented by diagrams where the number of 
external legs is (i.e., contributions to the total mean energy and the thermodynamical potential), 
2 (self-energy (Green's function)) and 4 (four-point vertex functions). 

The second step is an expansion of the regularized un-renormalized self-consistent self-energy 
around the self-consistent vacuum propagator which shows that the asymptotic behavior of the 
diagrams and sub-diagrams is ruled by their pure vacuum parts. 

Then an equation of motion for the temperature-dependent "infinite part" of the regularized 
self-energy is derived and it is shown that it can be renormalized by a temperature independent 
subtraction procedure. 

A. BPHZ-Scheme for the vacuum 



We first apply the BPHZ renormalization theorem p^, 22, pO[ for the vacuum. The only 
difference to the perturbative case is that we apply it to diagrams with self-consistent propagator 
lines. This is justified since Weinberg's power counting theorem is independent of the special form 
of propagators but only needs their asymptotic behavior stated above. 

We summarize the BPHZ scheme as follows. A sub-diagram 7 of a diagram T is defined as any 
set of lines and vertices contained in T which itself builds a proper vertex diagram: 7 C F. A sub- 
diagram 7 is called renormalization part if its superficial degree of divergence (or its dimension) 
is greater than or equal to 0. In our case this means it has at most four external legs. Two sub- 
diagrams 71 and 72 are called nested 71 C 72 if 71 is a sub-diagram of 72. If they have no line or 
vertex in common, 71 n 72 = 0, they are called disjoined. If they are neither nested nor disjoined 
they are called overlapping: 71 o 72. 

To any diagram T we denote the integrand following from the Feynman rules with I^. For a set 
of pairwise disjoined sub-diagrams 71, 72, . . . , 7n we write the integrand in terms of the integrands 
of the sub-diagrams /^. and the remainder of the integrand denoted by /r\{7i,...,7„}i usually called 



the reduced diagram: 



Ir = h 



r\{7l,-,7n} 






7j- 



(14) 



The original scheme by Bogoliubov and Parasiuk |21] defines recursively the integrand i?r of 
the renormalized diagram. If a diagram does not contain any renormalization part but is itself 
divergent it is called primitively divergent. In that case the renormalized integrand is defined by 
i?r = (1 ~ tY')Ir- Herein tr is the operator of the Taylor expansion with respect to the external 
momenta around up to the order of the dimension 5{'j) of the divergent diagram, which is in our 
case A — E: 



t'yl^iPl, 



,Pk) 



-5(7) , 

El, 



j=0 





E 



dU^{pi,... ,pk) 



fll,...,fMk>0 

^ll+...+^lk=j 
for 5(7) < 



dpi 



A*! 



• • dp't'' 



Pi'---Pk' 



P1-- 



=Pk=0 



(15) 



If the diagram is convergent the integrand is unchanged under renormalization. 

If the diagram is not only primitively divergent but contains divergent sub-diagrams the in- 
tegrand for the diagram with all subdivergences subtracted is called ^r and the renormalized 
integrand is defined by 



Rr 



Rr if S{T) < 

(1 - tr)^r if S{T) > 0. 



(16) 



From Weinberg's power counting theorem it follows that after this recursive procedure the integral 
over the internal momenta of Rr is finite. The definition of the counter terms by the Taylor operator 
t^ for any renormalization part 7 of the diagram shows that these are polynomials in the external 
momenta to the order 5(7) and thus can be interpreted as counter terms to the corresponding wave 
function normalization factors, masses and coupling constants in the original Lagrangian. 

Zimmermann solved Bogoliuobov's and Parasiuk's recursion with his forest formula. A forest 
is defined as any set of sub-diagrams (including the empty set and the whole diagram itself) which 
are pairwise non-overlapping. One can depict these sets by drawing boxes around the sub-diagrams 
and in a forest these boxes are not allowed to overlap but they can be nested. A forest is restricted 
if each of its boxes contains only renormalization parts. To each restricted forest 5" one associates 
again an integrand, namely 



^^ = n....(-*7)^r 



765^ 





_ir(4) 


1 1 


1 -ivw ^ 




1 _ir(4) . 




_ir(4) 


1 1 


_ir(4) 





(17) 



The diagram to the right shows an example case for a typical ladder diagram, which we shall 
consider in the following section. For this case (5(7) = for all sub-diagrams and the diagram itself 
such that for each box only the subdiagram value at vanishing external momenta is to be subtracted. 
The tilde over the product sign in (|l^) stands for the fact that in case of nested diagrams within 
the forest one has to apply the Taylor operators from the innermost to the outermost diagrams 
while for disjoined sub-diagrams the expressions are naturally independent of the order of Taylor 
operators, since then 

n 

Ir = -fr\{7i,...,7„} n ^Tfe- (^^) 

fc=i 

The forest formula then says that the integrand of the renormalized diagram is given by the sum 
over all restricted forests: 

It is understood that the empty set stands for the diagram itself, i.e., without any box around a 
sub-diagram. 

The described BPHZ-scheme chooses the renormalization point for the divergent diagrams at 
external momenta set to 0. It is clear that by another finite renormalization we can switch to 
any renormalization scheme appropriate for the application considered. In our case of (/)^-theory 
we chose the on-shell scheme. We have to define the coupling constant, the mass and the wave 
function normalization. This can be formulated in terms of the proper self-energy and the proper 
four-point vertex (the three-point vertex can be set identical to because of symmetry under field 
reflection <j) ^> —<j) without destroying the renormalizability of the theory, so that we do not have 
to consider terms linear or cubic in the fields within the Lagrangian): 

r(^'™^)(s,t,n = 0) = ^, S("^=)(/ = m^) = 0,9p2S(™^)(p2 = m^) = 0. (20) 

Here s, t, u are the usual Mandelstam variables for two particle scattering, p is the external 
momentum of the self-energy and rn^ is the renormalized mass of the particles. The first condition 
defines the coupling constant at vanishing momentum transfer for the two-particle scattering to 
be given by A, the second condition chooses m to be the physical mass of the particles, while the 
third condition ensures that the residuum of the propagator at p^ = m'^ is 1 and thus the on-shell 
wave function is normalized to 1 as it should be. 

B. The finite temperature self-energy at the regularized level 

In this section we like to isolate those vacuum subparts inherent in the pure temperature part 
of the self-energy which need to be renormalized at physical space-time dimension d = 4. For 
this purpose we assume a regularization scheme, e.g., dimensional regularization and extract those 
vacuum parts from the self-energy which diverge in the limit d — s- 4. For sake of clarification 
we mark all equations with an asterix which diverge in the limit d — > 4 and which need special 
renormalization treatment. All other equations are generally valid, even if all expressions are 
replaced by their renormalized quantities. 

In order to extract the divergent vacuum pieces we take the self-energy as functional of the 
self-consistent propagator and expand it around the vacuum value 

si2 = ssr^ + ss? + sg (21) 



^(matter) 
■"12 



Here E^™*^ is the vacuum (T = 0) self-energy. Its renormalization poses no particular problem 
and can be done according to standard rules. At the examples discussed in our second paper 
1 26 1 it is shown how to do this in practical cases for the numerical solutions to the self-consistent 
equations of motion. The second and third terms in (21) contain the in-matter or finite temperature 
components of the self-energy. Thereby 



-iS 



(0) 
12 



d(l'2') 



(5iS 



12 



5Gv 



T=0 



^(matter) 
"^1'2' 



cO 



ir(4) 



(22*) 



contains the parts of S linear in the matter (temperature) part (^('^attcr) q£ ^^iq full propagator 

iG 



•/^(vac) _|_ -/^(matter) 

+ 



(23) 



At this level it is important to recognize that all loops involving vertices from both sides of the 
real-time contour (cf. Fig. ||) are UV convergent due to the analytical properties of {— +} and 
{H — }-propagators (A16) to ( |A18D . Such loops contain at least one thermal weight factor n(po) 
which decays exponentially at large \p()\, since the 0-function parts completely drop out for large 
loop momenta. Therefore all mixed components like G '', or S"' have to be excluded from the 
subtraction scheme. Thus, the expansion point, G^^^'^' , in (22) is defined as the contour-diagonal 
part of the propagator in the vacuum limit (T — > 0), i.e., with vanishing G '"(™'^) = G"' '-^'^^ := 0. 
Likewise all vacuum structures like S'-™'^-*, and the four-point functions F^^'™'^^ and A'™'^^ defined 
below are "diagonal" in the real-time contour placement. 

The remaining self-energy piece T,^^' in ( ^I[ ) contains at least two Qi^^^^<^^) lines which therefore 
are never involved in any divergent loops due to the 2PI property of the $-functionaP . Thus, there 
are no hidden sub divergences in S'''\ and possible divergent vacuum sub-structures can directly 
be renormalized using the BPHZ rules given above. 

On the other hand the diagrams of S'^^^ deserve special attention, since there the single q('^^^^<^'^). 
line is involved in logarithmically divergent loops, if all vertices of F^g i'2' ^^^ placed on the same 
side of the contour. As mentioned the terms with mixed vertices are finite. 

The divergences result from the fact that the functional variation of S with respect to G at 
T = defines a vacuum vertex function 



-iF 



(4) 
12,1'2' 



5S 



12 



5Gii2' 



-2i 



52$ 



T=0 



0G120G1/2' 



(24) 



r=o 



with four external legs. 

Its diagonal part (all vertices on one contour side) defines F^2'™2' ^^ich is of divergence degree 0. 
Assuming (^l^attcr) q£ divergence degree —4 it is involved in a logarithmically divergent loop. Thus, 
this part of S^^^, called S^'^''^'^', accounts for all terms of divergence degree and consequently S^''^ 
is of divergence degree —2. 

In order to trace all subdivergences hidden in S'-''-*, the vacuum structure inherent in (^{matter) 
has to be resolved. For this purpose the diagonal parts of the full propagator (i.e., G and G^^) 



Any G'™''"'''''-line in E^''^ is either involved in loops with fiirther G'™''**"'^' -lines which are finite or it is attached 
to a pure vacuum piece. Due to the 2PI property of $ this vacuum piece has more than four external legs also 
leading to finite loops for this G^'^^'^^'^-line. 
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need further to be decomposed 



iG 



12 



iG 



(vac) 
12 



+ i j d(l'2') G^^^'^T.f4"^G^;f + iG 



(r) 

12 



|/^(mattcr) 



(25) 



+ 



_ip(4,vac) 



+ ^T5M^ , 



(26*) 



where here and below all vertices 1,2, 1' and 2' are placed on the same side of the contour, this 
way defining the divergent piece of H^^K The remaining part G^^' is of divergence degree —6. 

Since r^^'™^) is a pure vacuum four-point function, eqs. (21) and (|25|) also show S^'^''^'^) as a 
functional linear in cCmatter) 



.(0,div) 
"12 



|/r,/N ■p{4,vac) 



//r,//\ ^(vac)^(0,div)^{vac) ^(r) 



•^VJ^ ^ J ^ 12,1'2' \ / V / iJ^i/i// ^i//2" '-^2"2' "•" '-^1'2' 



.(r) 



and thus also linear in Gy2i- ^ simple iteration argument, starting with S 



(O.div) 
12 



(27*) 



0, shows that 



Eq. (|27*D is solved by the ansatz 



.^(0,div) 



/ d(l 2 ) Aj^2,l'2'^l'2' 

Jc 



0% 



jA(vac) 



(28) 



leading to the pure vacuum equation of motion 

. (vac) _ p(4,vac) . f j/nA ce^ -p(4,vac)^(vac)^(vac) . (vac) 
^'-12,1'2' ~ ^ 12,1'2' "•" '■ / Ul^O^dO; i ;^2 34 ^^35 ^^46 ^»-56.1'2' 

Jc 



(29*) 



for the vacuum four point function A^™'^^ The diagrammatic interpretation shows that this has 
the form of an inhomogeneous Bethe-Salpeter ladder equation 
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_jp(4,vac) 











(30*) 



From the construction it is clear that this is a very particular BS-equation, namely the one which 
complies with the self-consistent Dyson resummation scheme defined through (|ll|). Thus these 
ladders, which are of s-channel type (forward scattering), are implicitly contained in the self- 
consistent vacuum self-energy. Since Weinberg's power counting theorem shows that A'™'^^ is of 
divergence degree 0, as any four-point function, Eq. (^) implies that S''^' ^^-^ is indeed also of 
divergence degree as assumed above. Then G^^' is of divergence degree —6 and the loops in 
(pS), which close A'™'^) with a G^''^-line, contain no further divergences. Thus, once A'-™'^) is 
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renormalized, the self-energy is renormalized too, within a temperature independent subtraction 
scheme. 

Since $ is 2PI, the variational relation ( p4[ ) defines the BS-kernel ri2'™2' ^^ ^ proper skeleton 
diagram, i.e., it contains no self-energy insertions and cutting the diagram such that the pairs of 
space time points (12) and (1'2') are separated, cuts more than two lines. Thus the BS-kernel has 
the appropriate irreducibility properties for the resummation to the complete four-point function 
^(vac)^ again showing the virtue of the ^-functional formalism to avoid double counting. 

Prom (^) it is obvious that switching to the momentum space representation both r^^'™'^) and 
^(vac) g^j^g j^q|. ]2gg(jg(^ in their full momentum dependence but rather only as a function of the two 
momenta given by the Fourier transformation with respect to the space-time point pairs (12) and 
(1'2'). Due to (H) and through (p9*|) both, T^^'™'') and A(™'=), obey the symmetry relations 



p(4,vac) 
■*- 12,1'2' 
A (vac) 
^'-12,1'2' 



-p(4,vac) 
^ 1'2',12 
A (™c) 
^'-1'2M2 



or 



or 



r(4'™^)(p,g) = r(4=vac)(g^p)^ 

A(™c)(p^g) = A(™=)(g,p). 



(31) 



C. Renormalization of the vacuum Bethe-Salpeter equation 



In energy-momentum representation the regularized BS-equations ( p9*D and its equivalent "ad- 
joint" version become 



A(™^)(p,g) = r(4'™=)(p,g) + 



{2'kY 



r(^'™^)(p,0[G(™")(0]^A(™")(/,g) 



r(4'™^) (p, q)+\j (1^ A(-^) (p, /) [G(™^) (0]'r(4-™^) (l,q). 



(32*) 
(33*) 



The renormalization of the BS-equations ( |29*| ) is not straight forward. First the BS-kernel 
-p(4,vac)^^^ ^-j has to be renormalized following the BPHZ-rules outlined in sect. Ill, Represent- 
ing A(p, q) as the sum of ladder diagrams, this BS-kernel forms the rungs in each ladder. The 
complication arises from overlapping sub-divergences: in the general case each internal rung is 
part of various diverging sub-diagrams through the loops involving two or more rungs. However, 
two observations help to settle this renormalization issue: 

(i) The rungs given by the BS-kernel are 2PI with respect to cuts 

separating the top from the bottom extremities of the rung. 

This implies that there are no divergent sub-diagrams which 

cut into the inner structure of any rung, since this would involve 

a cut of more than two lines and the resulting sub-diagram 

would have more than four external lines. Thus divergent 

sub-diagrams always have complete rungs as sub-diagrams, as 

shown to the right. 

(ii) If one takes the difference of two A functions which differ only 
in one of the momenta, e.g., A(p, q) — A{p', q), all those counter 
terms cancel out which contain boxes which cut the outer p 
or p' lines, respectively, since for these counter terms the ar- 
gument p or p' are replaced by zero. Thus, the only boxes 
which are left are those which exclude the outermost rung at- 
tached to the p-lines. These boxes, however, just define the / i- i . . \ 
renormalized result of the sub-diagram complementary to this 
outermost rung, which is again a ladder diagram. 











_ir(4) 










1 




1 




_ir(4) 












_ir(4) 












_ir(4) 




1 




1 










-ir(4) 
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This permits to establish recursive relations for the renormalized expressions of the two possible 
differences 



A^'''''^ {p, q) - A^'™) (p', g)r(^'™=) {p, q) - r(4'™^) {p', q) 



+ i 



(2vr)^ 



[r(^'™'=)(p,/) - r(^'™'=)(p',/)][G(™'=)(/)]2A(''^'^)(/,g), 






(34) 



(35) 



+ 1 



(2vr)^ 



a('''='^)(p', 0[G^"^'^(0]^[r^^'™'H^ 9) - r^^'™") (/,(/')], 



where now r'^'™'^) stands for the renormalized Bethe-Salpeter kernel. Since the renormalized 
function T^^''''^^'{p,l) — r'^'™'^)(0, /) is of divergence degree less than zero the integrals are finite. 
This set of renormalized equations can be used to construct the renormalized A-function using 



a(''™)(o,o) 



(36) 



on the two real-time branches C^ due to our renormalization condition (^0|) . In a kind of sweep-up 
sweep-down scheme first the "half sided" A^^'^^^>{0,q) function can be constructed by solving (35) 
for p' = 0. Using this half sided function as the input for ( |3^ the full momentum dependence 
of A^''^") {p, q) can be obtained. This scheme fully complies with the BPHZ renormalization pre- 
scription. It has the remarkable feature that, although it is explicitly asymmetric in p and g, it 
constructs a completely symmetric renormalized four-point function which can be combined to the 
complete result 



A('--)(p,g)=r(4'™^)(p,g) 
dH 



+ 1 



+ i 



(2vr)4 
(2vr) 



[r(^'™=)(p, - r(^'™=)(o, /)][g(™^)(0]^a(™=)(/, q) 
^ A^'^""^ (0, /) [g(™^) (/)]2 [r(4'™^) {i,q)- r(^'™^) (/, 0)] . 



(37) 



For numerical applications it is important to realize that only the half sided A'''™^ (0, 1), and not the 
full momentum dependence A^'^'^^'{p,q), is explicitly needed. Since the half side A has essentially 
two-point function properties it can be numerically constructed using similar techniques as for 
self-energies. 

Indeed one can express the complete renormalized self-energy part linear in (^(matter) ^^ ^^^ form 



sW(p) = s(o)(p)-s(o)(o) + sW(o) 

Ml 

;^[T^^\p,i) - r(^'™=)(o,/)]G('^^"'''')(0 
+ / ^-^,^a('-™)(o,/)g('-)(/) 



{InY 




(38) 



-ir(4) 



_;p(4,vac) 



+ 



jj\(ren) 



-P 0- 



finite 
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with G^'"'""*^'') and G^'') from (H) and (|2|). Here r('*) is the full contour valued kernel (||) including 
mixed contour vertices, while r^^'™'^^ and A^''^'^^ are diagonal in the contour vertices. Due to the 
2PI-properties of T^^' the difference T^^'{l,p) — r(^'™'^'(/, 0) is of divergence degree less than 0. 
Therefore the first integral represented by the difference of the first two diagrams is finite, since 
^(matter) jg q£ divergence degree —4. This difference represents the most naive subtraction, which by 
itself, however, would be false, since it contains temperature dependent counter terms. The heart of 
the above derivation is that these false T-dependent counter terms are precisely compensated by the 
last term. The fact that the counter term structures never mix the two real-time contour branches 
also lifts the problem of pinch singularities which otherwise could arise due to the vanishing external 
momentum. 

This completes the proof that the self-energies can be renormalized with T-independent counter 
terms. 

D. Renormalization of the real-time F-functional 

In this section we derive the renormalized real-time generating functional for the in-matter 
equations of motion. Thus, we restrict the contour integrations to the real-time contour Cr and 
the corresponding traces to the real-time traces TriR. For functions in momentum representation 
the corresponding contour matrix algebra in the { — h} notation (cf. Appendix]^) is implied. 

For the renormalization procedure we use the ansatz (^) for the full propagator G together 
with the form (28) for the logarithmically divergent part of the self-energy, where the renormalized 



four-point function A resolves the sub divergences hidden in both, the propagator and the self- 
energy. For this purpose we decompose the generating functional F in its vacuum part, which is 
solely a functional of G^™'^^ and the in-matter part 

■p p(vac) r^(vac)"| _, p(matter) r/^(vac) /^(matter)] ^"^Q^ 

Thereby it is implied that the vacuum problem is already solved through its equation of motion 
resulting from the functional variation of F^™'^). Given G^™'^) the equations of motion in matter 
result from the functional variation of r(™'***°'') with respect to Cl^atter) ^ 

Compared to the two-point self-energy the F- and $-functionals have no external points and 
essentially result from the diagrams of the self-energies by closing the extremities. Therefore one 
has to explicitly expand the corresponding expressions up to second order in G(™atter) i^gfgj-g Qj^e 
comes to the situation where the remaining pieces are void of hidden sub divergences. Thus we 
write 



, , A(f)("^^c) 1 i:2ff)(vac) 

a„ - <T,(™'=) _L Tr-n= — /7(mattor) , ^ ^(matter) » ^ ^(matter) , ^(r) 

^" »^ + %G(--) ^2!^^ 5G(--)' "" (40) 

(J)(^^'^) _ T'y /-~i(mattcr)yi(vac) _, _/^(inattcr)p(4) ^(matter) _|_ (jj(i') 

Here we have used (0) and (^) for the vacuum parts defined through the variation of $ with 
respect to G. At the same time we introduced the real-time trace in momentum space 

Tr^A...B = I^Ail)---B{l) (41) 

and the functional tensor contraction for four-point functions with propagators: 

^^^^'^^^ = IW?I ^.G,{h)ri'Hw,l,)G,{k). (42) 
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In all expressions the functions are contour matrix functions which imply the corresponding contour 
matrix algebra and the contour trace, cf. Eqs. (A1S|) and (A20). 



Applying the arguments given for S^''^ also for <^'''' no G^™'**''°'''-line is involved in divergent 
loops such that after renormalization of possible vacuum sub-divergences the entire diagram is 
finite. Thus only terms with at most two G'^™^**°'')-lines need further care. 

Since ^ by itself is not an observable we directly step to the definition of the T-functional 
which relates to the thermodynamic potential in the equilibrium case. Exploiting the stationarity 
condition of T at the vacuum level, i.e., using the vacuum equations of motion for G'-™'^-*, all terms 
linear in c(™a.ttcr) ^^-q^ q^|^ g^j^^j -^g f^j^^^j fQj. ^i^q functional with all proper vacuum sub-divergences 
subtracted 

■ / °° //^(vac)y (matter)\fc \ 
p(matter)r^(matter)i __}_ rry [ /^(matter) v (matter) _ V^ I*-" ^ )__ \ 

^ V fc=2 ^ / (43*) 

-^TrKG('^^"'='')sW+'lJ;'\ 

Renormalized to zero at the vacuum level this expression is a functional of the in-matter part of 
the propagator. The remaining divergent parts of f (matter) g^j-jgg fj-om terms quadratic in c(inatter)^ 
i.e., 

p{matter,div) ^ i -Q.^ (^(0)5^(0) ^ l(-<(0)p(4,vac)(^(0) _^ l(-,(0)p(4,vac)(^(r) ^^^4*^ 

with 

QiO) _ (^(vac)j-j(0,div)^(vac) /^g^ 

where again in both relations above all expressions are contour diagonal. Both, S^"''^''^) and G^^' are 
linear in C'imattcr)^ "Qgjj^g the equations of motion for J]^''' '^^ and A^^^' one arrives at an expression 
for the divergent part of Yi^^^^'^'^i which only contains quantities which were already renormalized 
in the previous subsection 

-p(mattcr,div,rcn) _ 1 ( Qi^) J^i^cn) Q{r) _ ^(r)-p(4,vac)^(r)\ /^g-v 

Substituting this for the divergent part we obtain after some algebraic simplifications 

-p (matter, ren)r^(matter)-] _^tT;~/ ^(r)y (matter) _ ^^(vac)y(r)^(vac)y(r) 



fc=3 / 

+ mixed contour TrR-terms from ( [43*D , 



(47) 



where Tr^ includes only the contour diagonal parts. This expression, which now can be considered 
as a functional of (^(ma'tter)^ ^^ through Eqs. (p^) and ( p8| ) of G^"^', is void of any hidden subdiver- 
gences, since all matter or T-dependent parts of the propagator like Qi™^^^°^) or G^'^' are involved 
in convergent loops. 

Now it remains to be proven that this renormalization procedure for the T-functional is con- 
sistent with the renormalization of the self-energy given in the previous section. In other words: 
We like to show that the vanishing functional variation of p(rcn)j(^(mattGr)j complies with the Dyson 
equation of motion and the renormalized self-energy. Erom the BPHZ-formalism we expect this 
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to hold true, because it ensures that sub-divergences can be renormalized first and then the re- 
maining divergences which come into the game by closing the diagrams: The renormalized result 
is independent of the order of counter term subtractions. 

It is sufficient to show this for the contour diagonal parts of (0)- Writing the functional 
variation of r^"^*^"^ as 

i--p(niatter,ren) i--p(niattcr,ron) 
^ (matter,ren) ^ OL ^ (,) OL ^ (,) 

where S'''' is supposed to be a functional of G^^K Both terms in ( |4^ ) independently vanish. The 
first term drops by virtue of the Dyson equation which together with (25) ensures 



/ °° /'/^(vac)yi (matter) \fc \ , , 

S Y^ ^ L = (g^^) _ (-(vac)5.{r)^(vac)\ ^^.(mattcr)^ (-^g) 



\fe=3 
The second term shrinks to 



^j.(matter,ren) ^ l_ ^^^ (^^{r) ^Q{r) ^ f^{..c) ^{r) ^{..c) ^j.{0)^ ^^^W ^ q^ ^^q^ 



which indeed implies 



S(r)5G(™^"<='') 



s^'^:^ 



SW = 2i , ^^ , , (51) 

(^(^(matter) ' "^ ' 



compatible with the definition of S^'^^ It is important to note that through the functional variation 
( |4^ only convergent loops are opened, such that none of the counter terms is affected by this vari- 
ation. This explicitly demonstrates the consistency of the BPHZ-renormalization scheme for the 
self-consistent approximations: The operations of variation with respect to G and renormalization 
are commutative, i.e., one can construct the renormalized self-energy in two equivalent ways: The 
first uses the un-renormalized ^-functional and defines the renormalized self-energy by applying 
the BPHZ-renormalization theorem to its diagrams which are defined by opening any line of the 
un-renormalized ^-functional. In this way we have defined the renormalized self-energy in the pre- 
vious section. Subsequently we renormalized the F-functional by substituting the thereby defined 
renormalized functions for the divergent vacuum-sub-diagrams with 2 and 4 external legs. After 
subtracting the pure vacuum contribution this leads to the finite renormalized functional (^7|). 

As we have seen now, the second way to define the renormalized self-energy is to renormalize 
the F-functional first. Then the variation with respect to the renormalized propagator leads to the 
renormalized equations of motion and thus directly to the renormalized self-energy which shows 
the consistency of the local vacuum counter terms including all combinatorial factors for both, the 
generating functional and the equation of motion. 

We like to clarify that the equations derived in this section for the in-matter parts of the 
propagator are all valid also in the general non-equilibrium case of quantum field theory provided 
the density operator at time to amends a Wick decomposition. This is valid for statistical operators 



at initial time of the form [31 



R = — exp(- ^ OfcAfc) with Z = Tr exp(- ^ a^Afc), (52) 

k k 

where the A^ are one-particle operators. Our arguments for the renormalizability with Ofc inde- 
pendent (i.e., state independent) local counter terms should hold, since the statistical operator is 
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normalized TrR = 1. Thus, any in- matter part of the propagator leads to similar reductions of the 
degree of divergence as for the Bose-Einstein distribution functions used here. This ensures that 
the power counting arguments for the non-vacuum parts are still valid. Note in particular that 
the real-time functional T has only functional meaning, namely as a tool to derive the equations 
of motion, since its value at the physical solution vanishes. 

E. Renormalization of the thermodynamical potential 

Using the thermodynamic part, i.e., the vertical branch, of the contour given in Fig. |l| the T- 
functional provides a finite value which indeed relates to the thermodynamic potential Q. For the 
evaluation one uses the relationship between the Matsubara functions and the real-time functions 
given in Appendix ^, Eq. ( |A21 ) for the thermodynamic trace Trxh which takes due account of the 



thermodynamic weights in the partition sum 

= -TipV / j—^ sign(po) {n{po) + -jImhR{p), 

where h ^ and h^ are the Wightman functions and /i/j is the retarded function of h. Furthermore 

'^iPo) = 7Wi — IX 7 (54) 

exp(/3|po|) - 1 

is the thermal Bose-Einstein factor resulting from the summation over Matsubara frequencies 
expressed in terms of complex contour integrals cf. ( [A21 ). For this thermal contour Cxh closed 



diagrams as those of <& and F also attain a finite value. In this case the rule is, first to omit one 
of the momentum integrations, which in this way defines a two-point function. Its renormalized 
retarded value can be calculated according to the above used real-time contour rules. Subsequently 
one applies ( |53[ ) for the final integration. We show now that after renormalization of the two-point 
function we need only to subtract the overall vacuum divergence inherent in this final integral. 

To obtain this result we have to go back to the regularized expression for the un-renormalized 
effective potential @ which relates to the thermodynamical potential via 

j^('-^g)(r) = -Tr^'°s) (55) 

where 



i , / M2 \ i 



r5'™)(r)=TYTh 



-ln( -^^^ ] +1S™=G™" 
2 \ G™W 2 



p(,matter,regj ,rp\ _ •■ rp 



+ «&S?(T), 

(56) 



1 °° //^(vac) V (matter) \ A 

Ai(r)j^(mattcr) _ ;)^^{vac) j^(r)^(vac) j^(r) _ V^ V^ ^ ) 

9 2-^ h 



k=3 



Here the subscript Th specifies the quantities resulting from the thermal trace. The matter part 
results from form (|47|) . Subsequently one replaces all quantities by their renormalized ones (denoted 
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by a bar across the functions) and cancels the overah divergence by subtracting the T -^ +0 value 



f^(rcn)^^) 



=,(vac,rcg) 
Th 



r 



(reg) 
Th 



(+0) 



-T r 



(T)-r 



(vac, reg) 
Th 



(+0) + r^. 



(matter, reg) 



(T) 



(57) 



This procedure is legitimate as long as the new loops due to the final therniodynamical trace 
(|5^ do not induce new subdivergences, but only overall divergences. It is obvious that for the 
matter part the final trace loop involving the factor n + 1/2 is completely convergent, since all 
loops are regular once A is renormalized. For the two vacuum terms the component proportional 
to the factor 1/2 in the thermal trace (53) cancel out, such that all terms are proportional to n{po) 
which cuts off the po integration, while the momentum integrations are also limited due to the 
vacuum thresholds: the imaginary parts of vacuum functions are zero for p^ < m^, where m is the 
mass of the stable vacuum particle. Thus also these final loop integrals are finite, defining a finite 
thermodynamical potential. 

The vacuum part essentially determines the kinetic energy part of $7 as can be seen from the most 
simple example of an ideal gas. Here of course all self-energies and $ are vanishing, the retarded 



propagator at finite temperature is -Dr(p) = [p 



m' 



' -|-i?7cr(po)] ^-nd the renormalization is done 



by subtracting the pure vacuum part. Thus, the free thermodynamical potential becomes 



j^(id. gas) ^ _Y 



Si 



■n(/o)^e(/2-m2), 



which can be brought to a more familiar form by an integration by parts 



J^(id. gas) ^ _py ^y I 



d^ 



(2vr)= 



In 



1 — exp 



-HP 



+ m' 



(58) 



(59) 



IV. CONCLUSIONS AND OUTLOOK 



For the example of (/> -theory we have shown that self-consistent Dyson resummations based on 
a ^-derivable scheme can be renormalized with local counter terms defined on the self-consistently 
determined vacuum level. This result was obtained with help of Weinberg's power counting theorem 
and using the BPHZ-renormalization scheme with the usual modifications for finite temperature 
diagram rules, which can be summarized in the simple rule that the "contraction boxes" defining 
the counter terms have to exclude sub-diagrams which contain any temperature line. 

The hidden subdivergence structure of the self-consistent scheme has been resolved. This leads 
to a Bethe-Salpeter equation for the vacuum four-point function compatible with the chosen $- 
approximation, which we have renormalized. The method is free of pinch singularities. Closed 
equations could be formulated which resum the non-perturbative structure of both, the equations of 
motion, i.e., the self-energies and also the non-perturbative counter-term structure. The complexity 
of these equations is comparable to standard Dyson resummation schemes and therefore in principle 
does not imply new techniques. First numerical applications, which include the construction of 
the BS-kernel, the solution of the half-sided four-point function and thus the renormalized self- 
energies up to the self-consistent sunset self-energy, are presented in a second paper |26|. The 



renormalization of the generating functional F, c.f. sect. HID, shows that the derivation and thus 
the renormalized in-matter equations of motion equally apply to the general non-equilibrium case. 
This also proves that there is no arbitrariness in studying the in-medium modifications of 
model parameters like the mass and the coupling constants within this class of approximation 
schemes: It is sufficient to adjust them in the vacuum, for instance by fitting them to scattering 
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data, in order to predict without ambiguity how they change in the dense and hot medium: The 
in-medium m,odifications are ruled completely by the model and its vacuum parameters alone, no 
further assumptions need to be made. 

Although demonstrated for the i;^^-theory, the method is in principle general, since the derivation 
only relies on the analytic and asymptotic form of the propagators. In particular the renormaliza- 
tion of hidden overlapping divergences in the logarithmically divergent Bethe-Salpeter equations 
is general. Still, there is a number of restrictions of the self consistent Dyson resummation within 
the ^-derivable scheme, which concerns global and local symmetries and the corresponding con- 
servation laws and Ward-Takahashi identities. 

The ^-functional formalism only ensures the conservation laws for the expectation values of 
charges associated with the symmetry by Noether's theorem. However, in general the Ward- 
Takahashi identities are violated for the self-energy and higher vertex functions. Heuristically the 
problem can be traced back to the violation of crossing symmetry by the self-consistent scheme: 
Our derivation shows that the self-consistent solution of the self-energy involves Bethe-Salpeter 
ladder resummations of the four-point function, but only in the s-chanel. The crossing symmetric 
t- and u- channel contributions to the four-point function are not included. 

The symmetry properties of the r[99, G]-functional were already investigated by us with help 



of the here applied path-integral method [^, 32]. We show that it is always possible to define 



a non-perturbative approximation to the effective action reff[</?] which respects linearly realized 
symmetries of the classical action provided the symmetry is not anomalously broken. The self- 
energy and higher vertex functions defined from this improved approximation action formalism 
then fulfill the Ward-Takahashi identities of the underlying symmetry. As a result the effective 
action Peg enforces that additional t- and u-channel Bethe-Salpeter resummations are needed to 
restore the crossing symmetry together with the Ward-Takahashi identities for the self-energy and 
the vertex functions. However, these vertex functions are not self-consistently calculated and thus 
some problems remain also within this approximation: For instance in the case of the linear sigma- 
model the 0(N)-symmetry is restored for the vertex- functions and the Goldstone- modes become 
massless. Yet, the phase transition from the Nambu-Goldstone phase at low temperatures to the 
Wigner-Weyl phase at high temperatures results to be of 1^* order rather than 2°*^ order Q. 

In the case of a local gauge symmetry the problems become even more intricate: Self-consistent 
schemes beyond the classical field level for the gauge fields generally violate local gauge symmetries 
for the same reasons as for global symmetries. However, this immediately causes the excitation of 
spurious modes of the gauge fields which leads to violation of the unitarity of the S-matrix, the 
positive definiteness of the statistical operator and the causality structure of Green's functions. 
Nevertheless a gauge invariant effective action Fcfr within a background-field approach can be 



formulated which provides gauge covariant polarization functions [17|. 

From a practical point of view the problem remains to calculate the self-consistent propagators 
needed for the symmetry-restoring Bethe-Salpeter resummation, which presently can only be solved 



in simple cases (RPA bubble resummation). In |17] we have presented a workaround in terms of a 



suitably chosen projection method onto the physical (transverse) degrees of freedom of the gauge- 
field polarization tensor. This procedure, of course, does not lead to a full restoration of local 
gauge theory but to causal Green's functions and current conservation within the self-energies of 
matter-fields. Alternative methods are restricted to the approximate solution of the self-consistent 
equations of motion, e.g., in the sense of a Hard thermal loop approximation |19] or to a systematic 
expansion in terms of the coupling constant or h [^ . 

The proof of the renormalizability of ^-derivable approximations opens a broad range of per- 
spectives for effective field theory model applications describing the non-perturbative in-medium 
properties of particles in dense or finite-temperature matter with model parameters fixed at the vac- 
uum level. Further applications point towards the appropriate renormalization of non-equilibrium 
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transport equations [^], where in particular the drift terms, which determine the equation of state, 
involve the real part of the self-energies which generally need renormalization. 
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APPENDIX A: ANALYTICAL PROPERTIES OF GREEN'S FUNCTIONS 

In this appendix we summarize briefly the analytic properties of Green's functions of neutral 
bosons needed in the main part of the paper. This is most easily done by switching to the operator 
formalism in the Heisenberg picture. By definition we have for a hermitian scalar field operator 

iGix) = I Tr exp(-/3H)0(x)(/)(O) := (Tc0(x)0(O))^ . (Al) 

For x^ on the vertical part, i.e., x^ = —it with < r < /3 we obtain the Matsubara Green's 
function 

GM{r,x) = G+-{-iT,x), (A2) 

where one has to understand the analytic continuation of the real-time Wightman function G"' 



on the right hand side. It is important to keep in mind that (A2) is only valid when the first time 
argument x^ in ( [Aip is on the vertical part of the contour while the second one is at x^ = 0. If 
both arguments of the fields are located on the vertical part according to ( [Al| ) one has to use the 
time ordering along the imaginary time axis. 



Since the order of the operators under the trace in ( Al) can be changed cyclically the real-time 
Wightman functions G"* and G ^ are related through 

G+-(x°-i/J,f) = G-+(xO,f), (A3) 

where x^ is a real-time argument on the contour and on the right hand side one has to understand 
the analytic continuation. One should keep in mind that only the Wightman functions are ana- 
lytically continueable, not the time ordered or anti-time ordered of the upper or lower real-time 
branch, since they contain step functions from the time ordering operator Tq. The same holds true 
for the imaginary-time ordered Matsubara Green's function. 

As in vacuum quantum field theory due to translation invariance it is customary to use the 
energy-momentum representation of Green's functions. For the real-time propagators we have the 
usual description 

G'^ix) = / (^ exp(-ipx)G^^-(p). (A4) 

We write down the formalism for arbitrary space-time dimensions since all considerations do not 
depend on it and we need it to obtain well defined non-renormalized quantities in the sense of 



dimensional regularization. The periodic boundary condition (A3) translates into the Kubo-Martin- 



Schwinger- condition (KMS) for the Fourier transformed Green's functions: 

G-+{p)=expi-Ppo)G+~ip). (A5) 
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The Matsubara Green's function is only defined for imaginary times — ir with < r < /?. Thus 
the momentum representation with respect to the time component is a Fourier series with period 
P according to ( [A3[ ) rather than a Fourier integral: 



Gm{x) 



1 



E 



^d-i 



{2ttY 



j-exp(-ij5x)GM(ipo,R)lpO=<^„ with 



UJ.n 



27r 



-n. 



(A6) 



Herein x*^ = — ir with < r < /5. Using the inverse Fourier transform, Eq. ([A4|), and the 

with /3(p) = i[G+"(p) - G"+(p)]. 



KMS-condition (A5) we find the spectral representation 



GM(iw„,,R) 



27r po - it^r^ 



(A7) 



This shows that the Matsubara propagator is completely represented by the real-time function p. 
With help of this we define the analytically continued propagator by 



Gc{k) 



dpop{po,k) 
27r po - ^0 



It can contain singularities only on the real axis. Using the Fourier transformation (^ 
limits to the real axis from above and below we obtain: 



GciPo ±iv,P) = Gr/a{p), Po G K 
with the retarded and advanced Green's functions 

G,j/A(x) = Tie(±t)([0(x),0(O)]_>^. 
From this we find immediately 

Gr{x) = G\{-x) => Gr{p) = G*Ap), p{p) = -2ImGK(p) 



-p{-p) 



and from (A8) and the analyticity of Gc{p) in the upper complex po-plane it follows that 

(y{pQ)p{p) > 0. 

For later use we note the momentum space properties 

G~~ + G++ = G+~ + G~+, Gr = G"~ - G~+, Ga = G" - G+-, Gm = -iGc(iw„) 



(A8) 
for the 

(A9) 
(AlO) 
(All) 
(A12) 
(A13) 



which follow immediately from ([Al]), (|A7|), and ( [AlO ). 

We also make use of the expressions for the real-time Green's functions in terms of the retarded 
Green's function, which follow immediately from (A7, A11 , |A18 ): 



iG--(p) = iGR{p) + [e(-po) + n{po)]p{\po\,p), 
iG++(p) = [@{po) + n{poM\polp) - iGRip), 
iG-+(p) = [e(-po) + n{po)]p{\po\,p), 
iG+~(p) = [&{po) + n{po)]p{\po\,P)- 



and the Bose-Einstein distribution defined as 

n{po) = 



exp(/3|po|) - 1' 



(A14) 
(A15) 
(A16) 
(A17) 



(A18) 
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All relations given above for the Green's functions G directly apply to any two-point function given 
by local field operators ih{x,y) = (7c^H(x)H(y)), e.g., the self-energy. 

Real-time contour integrations and traces of translationally invariant two-point functions 



Cix,y)= dzA{x,z)B{z,y), x,y,zeC^ 



T¥rC= / dxC{x,x) 
Jck 



(A19) 



transcribe to 



C{pr =Y,Mpt<ykiB{pf, a = diag(l,-l), i,j,fc,? £{-,+} 



kl 



ij 



(A20) 



in contour momentum-space representation. 

We close this appendix by citing the formula for summation over the Matsubara frequencies 



needed when calculating quantities related to the vertical branch of the contour |24]. In this paper 
we use this to calculate the thermodynamical potential 



1 



Tr^^ h{p) := PV- J2 



j^:^^h{iu;^,p) 



^dp 



(27r)"' 



2 + "-(^o] 



-f3V 



{h[po + ir/o-(po)] - h[po - ir]a{po)]} 



(A21) 



In ( A21| ) we have assumed that the function h is analytic below and above the real axis and that 
the trace exists. Usually this is only the case for the regularized or the renormalized functional 
traces. It is also clear that due to the exponential damping from the Bose-Einstein distribution 
( A18| ) this part of the integral has a superficial degree of divergence reduced by 1 compared to the 
first part which is not damped by an n- factor. 

Eq. ( A21) shows that the thermodynamical potential can be calculated from real-time quantities 
since the analytic continuation of the Matsubara Green's function needed on the right hand side 
is unique and can be obtained from the retarded Green's function as well c.f. ( |A9D . For a more 
detailed analysis of the analytic properties see also |25| and for the general case of Wigner function 
representations in the non-equilibrium context ||ll||. 
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